%% function to provide source term for
% dy/dt + Beta*exp(y)= 0 

function [dc] = source_stiff(t,c_old)
beta = 2;
t=t+1;
t=t-1;
dc = - beta * exp(c_old);

return
end
